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Boundary Integral Equations for the Transmission 
Eigenvalue Problem for Maxwell’s Equations 

Fioralba Cakoni* Houssem Haddar^ and Shixu Meng^ 


Abstract 

In this paper we consider the transmission eigenvalue problem for Maxwell’s 
equations corresponding to non-magnetic inhomogeneities with contrast in electric 
permittivity that changes sign inside its support. We formulate the transmission 
eigenvalue problem as an equivalent homogeneous system of boundary integral equa¬ 
tion, and assuming that the contrast is constant near the boundary of the support 
of the inhomogeneity, we prove that the operator associated with this system is 
Fredholm of index zero and depends analytically on the wave number. Then we 
show the existence of wave numbers that are not transmission eigenvalues which by 
an application of the analytic Fredholm theory implies that the set of transmission 
eigenvalues is discrete with positive infinity as the only accumulation point. 

Keywords: The transmission eigenvalue problem, inverse scattering, boundary integral 
equations. Maxwell’s equations. 


1 Introduction 

The transmission eigenvalue problem is related to the scattering problem for an inhomo¬ 
geneous media. In the current paper the underlying scattering problem is the scattering 
of electromagnetic waves by a (possibly anisotropic) non-magnetic material of bounded 
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support D situated in homogenous background, which in terms of the electric held reads: 


curl curl E^ — k‘^W = 0 

in 

R^\D 

(1) 

curl curl E — k'^N'E = 0 

in 

D 

(2) 

1/ X E = 1/ X E^ UM X E* 

on 

dD 

(3) 

u X curl E = 1 / X curl E'^ -|- z/ x curl E* 

on 

dD 

(4) 

lim (curl E® x x — ikrW) = 0 



(5) 


r^oo 


where E* is the incident electric held, E® is the scattered electric held and N{x) = 
€(x^ (j(x^ 

- \-i - is the matrix index of refraction, k = u^eoHo is the wave number cor- 

Cq coeo 

responding to the background and the frequency cu and the Silver-Miiller radiation condi¬ 
tion is satished uniformly with respect to a; = x/r, r = |x|. The diherence iV — in the 
following, is refereed to as the contrast in the media. In scattering theory, transmission 
eigenvalues can be seen as the extension of the notion of resonant frequencies for impen¬ 
etrable objects to the case of penetrable media. The transmission eigenvalue problem is 
related to non-scattering incident helds. Indeed, if E* is such that E'^ = 0 then E|/) and 
Eq = E®|£) satisfy the following homogenous problem 


curl curl E — k'^N'E = 0 in 
curl curl Eq —/c^Eq = 0 in 
z/ X E = 1 / X Eq on 

u X curlE = u X curlEo on 


D 

(6) 

D 

(7) 

dD 

(8) 

dD 

(9) 


which is referred to as the transmission eigenvalue problem. Conversely, if ([^-(|^ has 
a nontrivial solution E and Eq and Eq can be extended outside D as a solution to 
curl curl Eq — fc^Eo = 0, then if this extended Eq is considered as the incident held the 
corresponding scattered held is E^ = 0. 


The transmission eigenvalue problem is a nonlinear and non-selfadjoint eigenvalue problem 
that is not covered by the standard theory of eigenvalue problems for elliptic equations. 
For a long time research on the transmission eigenvalue problem mainly focussed on show¬ 
ing that transmission eigenvalues form at most a discrete set and we refer the reader to 
the survey paper for the state of the art on this question up to 2010. From a prac¬ 
tical point of view the question of discreteness was important to answer, since sampling 
methods for reconstructing the support of an inhomogeneous medium H. i fail if the 
interrogating frequency corresponds to a transmission eigenvalue. On the other hand, due 
to the non-selfadjointness of the transmission eigenvalue problem, the existence of trans¬ 
mission eigenvalues for non-spherically stratihed media remained open for more than 20 
years until Sylvester and Paivarinta 27 showed the existence of at least one transmission 
eigenvalue provided that the contrast in the medium is large enough. A full answer on 
the existence of transmission eigenvalues was given by Cakoni, Gintides and Haddar 
where the existence of an inhnite set of transmission eigenvalue was proven only under the 
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assumption that the contrast in the medium does not change sign and is bounded away 
from zero (see also II H. [I4] and for Maxwell’s equation). Since the appearance of 
these papers there has been an explosion of interest in the transmission eigenvalue prob¬ 
lem and the papers in the Special Issue of Inverse Problems on Transmission Eigenvalues, 
Volume 29, Number 10, October 2013, are representative of the myriad directions that 
this research has taken. 

The discreteness and existence of transmission eigenvalues is very well understood under 
the assumption that the contrast does not change sign in all of D. Recently, for the scalar 
Helmholtz type equation, several papers have appeared that address both the question 
of discreteness and existence of transmission eigenvalue assuming that the contrast is 
of one sign only in a neighborhood of the inhomogeneity’s boundary dD, 0. H. 

181 and 


23 , 24 


29 


15 


The picture is not the same for the transmission eigenvalue 
problem for the Maxwell’s equation. The only result in this direction is the proof of 


discreteness of transmission eigenvalues in 10 for magnetic materials, i.e. when there is 
contrast in both the electric prematurity and magnetic permeability. The T-coercivity 
approach used in does not apply to our problem which mathematically has 

a different structure form the case of magnetic materials and this paper is dedicated 
to study the discreteness of transmission eigenvalues for the considered problem under 
weaker assumptions of N — I. Before specifying our assumptions and approach let us 
rigorously formulate our transmission eigenvalue problem. 

Formulation of the Problem; Let H G be a bounded open and connected region 
with C^-smooth boundary dD := T (we call it T for convenience of notation as will be 
seen later) and let u denotes the outward unit normal vector on T. In general we consider 
a 3 X 3 matrix-valued function N with L°°{D) entries such that ^ ■ Re(iV)^ > a > 0 and 
^ ■ Im(iV)^ > 0 in H for every ^ G C^, |^| = 1. The transmission eigenvalue problem can 


be formulated as Ending E, Eq G L^(Z1), E — Eq G Ho(curl^, D) that satisfy 

curl curl E — fc^iVE = 0 in D 

(10) 

curl curl Eq — fc^Eo = 0 

in D 

(11) 

z/ X E = Z/ X Eq 

on T 

(12) 

z/ X curl E = z/ X curl Eq 

on T 

(13) 

where 

lV{D) := 

= 1 . 2 , 3 }. 



H(curl^,Zi)) := {u : u G L^(Zi)), curl u G L^(Zi)) and curlcurlu G L^(Zi))} , 
Ho(curl D) := |u ; u G H(curl^, D), = 0 and y^curl u = 0 on T} . 


Definition 

L2(T)), e- 


1.1 Values o//c G C for which the (10)-(13) has a nontrivial solution E, Eq G 
Eq G Ho(curl^,Zl) are called transmission eigenvalues. 
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It is well-known 16 that, if Re(iV — I) has one sign in D the transmission eigenvalues 


form at most a discrete set with -|-oo as the only possible accumulation point, and if in 
addition Im(iV) = 0, there exists an inhnite set of real transmission eigenvalues. Our 
main concern is to understand the structure of the transmission eigenvalue problem in 
the case when Re(iV — I) changes sign inside D. More specihcally in this case we show 
that the transmission eigenvalues form at most a discrete set using an equivalent inte¬ 
gral equation formulation of the transmission eigenvalue problem following the boundary 
integral equations approach developed in 15 . The assumption on the real part of the 


contract N — I that we need in our analysis will become more precise later in the paper, 
but roughly speaking in our approach we allow for Re(A^ — I) to change sign in a compact 
subset of To this end, in the next section we consider the simplest case when the 
electric permittivity is constant, i.e. N = nl with positive n 7 ^ 1, for which we develop 
and analyze an equivalent system of integral equations formulation of the corresponding 
transmission eigenvalue problem. This system of integral equations will then be a building 
block to study the more general case of the electric permittivity N. We note that the 
extension to Maxwell’s equations of the approach in 15 is not a trivial task due to the 


more peculiar mapping properties of the electromagnetic boundary integral operators as 
it will become clear in the paper. 


2 Boundary Integral Equations for Constant Electric 
Permittivity 

Let n > 0 be a constant such that n ^ 1 and consider the problem of hnding E, Eq G 


L^(Zl), E — Eq G Ho(curl^, D) that satisfy 

curl curl E —/c^nE = 0 in D (14) 

curl curl Eq — fc^Eo = 0 in D (15) 

z/xE = i/xEo on T (16) 

u X (curlE) = u X (curlEo) on T (17) 


In the following we set ki := k^/n. Before formulating the transmission eigenvalue problem 
as an equivalent system of boundary integral equations, we recall several integral operators 
and study their mapping properties. To this end, let us dehne the Hilbert spaces of 
tangential helds dehned on T: 

m’'’2(div,r) := {u G Hr(r),divru G 
H*i’"2(curl,r) := {u G (r),curlru G H"’^(r)} 

endowed with the respective natural norms, where curlr and divr are the surface curl 
and divergence operator, respectively, and for later use Vr denotes the tangential gradient 
operator. (Note that the boldface indicate vector spaces of vector helds, whereas non-bold 
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face indicate vector spaces of scalar fields.) If 7 r u = i/ x (u x z/) denotes the tangential 
trace of a vector field u on the boundary F, we dehne the boundary integral operators: 


Tfc(u) := 


= ^7r(fc^/ <hfc(-,y)u(y) dsy + Vr j y)divru(y) , (18) 


and 


where 


Kfc(u) := 7 r ( curl / <l)fc(-, ?/)u(y) dsj 


^ ^ik\x-y\ 


(19) 




dvr \x — y\ 

is the fundament solution of the Helmholtz equation Au + k‘^u = 0. Referring to 15 


and 25 for the mapping properties of the single layer potential 


X(v>) := / 


( 20 ) 


with scalar densities (^, we have that the boundary integral operator 


Sfc(u) = / $fc(-,y)u(y)ds 


( 21 ) 


acting on vector helds u is bounded from H~ 2 +®(r) to H 2 ’''®(r) for —1 < s < 1 and hence 

Tfc : “5(div, F) -)■ -i(curl, F) 

Kfc ; H^^’^5(div, F) —>■ (curl, F) 


are bounded linear operators. Now from the Stratton-Chu formula 11 we have that 


Eo(x) = curl y (Eo X z/)(y)$fc(x,y)dsy + J {cm\Eo x u){y)^ki^,y)dsy 
+ j divr(curlEo X z/)(y)<l)fc(x,y)dsy for xeD 

with similar expression for E where k is replaced by fci := k^/n and hence we have 
the integral expression for E — Eq. Note by taking the difference E — Eq we have the 
corresponding kernel <l>fcj(x, y) — <l>fc(x, y) is a smooth function of x, y, and approaching 
the boundary F and noting E x v = EqX v and curl E x v = curl Eq x z/ we have 

7 r(E — Eq) = (Kfc — KfcJ(Eo x z/) + (curlEg x z/), 

7 rcurl(E-Eo) = (K^ - KfeJ(curl Eq x z/) + (A;Tfc - (Eq x z/). 
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From the boundary conditions (16) and we have 7r(E—E q) = 0 and yrcurl (E—Eq) = 

0, i.e. 


Kfc(Eo X I/) + iTA.(curlEo X z/) - Kfcj(E x z/) - ^Tfci(curlE x z/) = 0, (22) 
Kfc(curlEo X z/) + A;Ta;(Eo x z/) — Kfcj(curlE x z/) — /ciTfc^(E x z/) = 0. (23) 


Introducing M = Exz/ = EoXz/ and J = curlE x u = curlEo x z/, we arrive at the 
following homogeneous system of boundary integral equations 


f kiTk^ — kTk — Kfc \ / 0\ 

for the unknowns M and J. Let us dehne 


(24) 


f fciTfcj — kTk Kfci — Kfc \ _ f — kTk \ 


Note that while the operator is a smoothing pseudo-differential operator of 

order 2 (see e.g. 15 and (^), the operators in the main diagonal have a mixed structure. 


Indeed, from the expressions 


kiT,, - kTk = (klSk, - k^Sk) + Vr o {Sk, - Sk) o divr (26) 

— -Tfc = {Ski ~ Sfc) + Vr o ° 


where S and S are dehned by (20) and (21) respectively, we can see that these operators 
have different behavior component-wise. Hence a more delicate analysis is called for to 
hnd the correct function spaces for M, J and their dual spaces in order to analyze the 
mapping properties of the operator L(/c). 


Lemma 2.1 The dual space of ^ 2 ’ 2 (div, r)zs H-ii(curl,r). Foru* G H 2’2(curl,r) 
and u G H's “s (div, F), (u*,u) is understood by duality with respect to L^(F) as a pivot 
space. 


Proof. For any tangential helds u G H 2> 2 (div, F) and u* G H 2 >2 (curl, F), we con¬ 
sider the corresponding Helmholtz orthogonal decomposition 


u = curi Y(1 + VrP, = curl r<?* + Vrp*. 


26 


Since divru = divrVrP = ^rP ^ H 2 (F) we have by eigensystem expansion (e.g. 
that Vrp G H 2 (F). Similarly, from the fact that curlpu* G H 2 (F) we obtain that 
curlrg* G Hi(F) . Now 


(u*,u) = (^curirg* -h Vrp*, curlrg + VrP^ 
= curlrg*, curlrg^ -h (VrP, Vrp‘) 
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3 3 

Hence the right hand side is well dehned in the sense of duality of H 2 (r)-H~ 2 (r) and 
H2(r)-H-2(r), and thus H 2-2 (curl, F) is in the dual space of H 2 > 2 (div, F). 

Furthermore, if u* = curlrO'* + Vrp* is in the dual space of H~i’“5(div, F), then (u*, •) is 

3 1 -^ 

continuous and linear on H~ 2>“2 (div, F). Then for u = curlrQ' 

u) = (^curlr<?*, curfrO'^ • 

Notice curi yQ is only in H~5(F), therefore by eigensystem analysis curlrO'* ^ Hi(F) and 
curlrcurfrg* G H^(T), i.e. curlpu* G H^(T). Now for u = VrP where VrP G H5(F) 

(u\u) = (Vrp*, Vrp) • 

Then VrP* G H^5(F). Therefore u* G H~5’5(curl,F). Now we have proved the lemma. 

□ 


In the following the spaces H 


other in the duality dehned in Lemma 2.1 


2 (div, F) and H 2-2 (curl, F) are considered dual to each 
In the next lemma we establish some mapping 


properties of the operator L(fc) given by (25). 


Lemma 2.2 For a fixed k, the linear operator 

L{k) : H7^(F) X H'i’-5(div,F) ^ H|(F) x (curl, F) 

is bounded. Moreover, the family of operators L(fc) depends analytically on k E C\M_. 

Proof. Let E,Eo G L^(Zl), E — Eg G Ho(curl^,Zl) be a solution to the transmission 
eigenvalue problem (14)-((T7|). Hence 


M = E X z/ G Hi ^(F), J = curlE X 1 / G Hi "(F). 

Noting that divr(curlE x u) = curl rcurl E = curl^E ■ z/|r, we have that divpJ G Hi " (F) 

and therefore (M, J) G Hi "(F) x H~i’“3(div, F). It is known from 15 that S^, — S^, 

Kfci - Kk are smoothing operators of order 1, 3 and 2 respectively. Then using ( [2^ we 
have that the following operators are bounded 


— kTk 

_ A'T' 

fci fc 


HX(r)^Hf(r) 

HF(r)^Hf(r) 

H-t-5(div,r) -» Hd(r) 


Moreover 


curl 1 


-Kfc)M + (^Tfc,-iTfc)J 


= curlr(Kfc, - Kfc)M + curlr(Sfci - 8^)3 G HKF), 
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and hence 


(fciTfc, - fcTfc)M + (Kfc, - Kfc)J G H|(r), 

(Kfc, - Kfc)M + J G H-5’i(cnrl,r), 

Hence L(fc) is bounded. Note that since every component of L(fc) is analytic on C\]R_, 
then L(fc) is analytic on C\M_ (recall that ki = ky/n). □ 

We need the following lemma to show the equivalence between the transmission eigenvalue 


problem and the system of integral equations (24). 


Lemma 2.3 Let Vt he any bounded open region in and denote V(curl^,f2) := {u : u G 
L^(fl), curl^u G L^(f2)}. For (p G ^(F), -ip G H~i “5(div, F), we define 


and 


Mi((p)(x) := curl J <hfc(x,y)(p(y)dsj^, x G M^\F, 


M2(^/’)(y) := / ^k{:>^,y)^p{y)dsy, a: G ]R^\F. 


Then Ml is continuous from ^(F) to V(curP, and M 2 is continuous fromii~ 2 ’^ 2 (^div,T) 
to V(curl^, D^) where D~ = D and D~^ = Bji\D with a sufficient large ball containing 
the closure ofD. Furthermore the following jump relations hold 


[7iMi(v7)] = (f 

in 

KHT), 

(27) 

[7tCurlMi((p)] = 0 

in 

HGff), 

(28) 

[7tCurlM2(V')] = 

in 

Hd(r), 

(29) 

[divr7tCurlM2(V’)] = divr'^ 

in 

H-i(n 

(30) 


Proof. Let us denote by < •, • > the (F)-Hi ^(F) or H 2 {Y)-H 2 (L) duality product. 
_ 1 

Since (p G ^(F), then from the classical results for single layer potentials 


|Mi((y9)||L2(D±) < C 


^fc(x,y)v7(y)dsj 


< c 


Hi(D±) 


H, 2 (r) 


and since curl^Mi((p) — fc^Mi((p) = 0 in then 

||curl^Mi((p)||L2p±) = |/c2|||Mi((p)||l2(d±) < c 




3 1 

where c is some constant depending on k. For fir G H~ 2'“2 (div, F), we have from 


15 


||M2(^/>)||l2(D±) < cWfiW 3 

/ 









Notice that 


curl^M2(^/^)(x) = / $fc(x,y)?/>(y)rfsy + V / divr^Piy)^ki-,y)dsy 


and dlyrt/^ ^ H 2 (r), hence we have from 


||curl^M 2 (^/')||L 2 (n)±) < c 


15 


3 + ||divr'?/’ll 1 , . 

H72(r) " 


This proves the continnity property of Mi and M 2 . To prove the jnmp relations, we will 
use a density argument. Let 


= curl / ^ki^,y)ip{y)dsy in D 


± 


We dehne the tangential component by duality. For a E (T), ||a|| 1 =1, there 

(r) 

exists G and w+ compactly supported in Bji such that 7 iCurl w = 0, 7 tw = a 

and ||w||h2(_d±) < c||q!|| 1 (see 16 ). Moreover, 

Hf (r) 


< a,7tu^ >= ± / (u^ ■ curKw^ — ■ curmu^)(ix. 

Jd± 


Then 


I < a!,7iu=^ > I < (||u||l2(d±) + ||curl^u||L2(D±))||w||H2(D±) 

< Ci(||u||l2(i 5±) + ||curl^u||L 2 (D±)) 

< C2||(p|| _1 

H, ^T) 

where ci and C2 are independent from u, therefore ||7iu=*"|| _i < C2||(p|| _i . Choosing 

Hj 2(r) ” 2.T.N 

2 ^ 

ipn E H~2 “2(div, T) such that —>■ 93 in ^(T) yields 

||7iU^ - 7tU^IL-1 < c||v? - -1 ^ 0, 


H, cn 






since [ 7 tU„l = (pn for (p„ G H 2 ’ 2 (div, T) (see 26 ). Letting n —)■ cx) yields [ 7 iul = (p in 

—- r 1 *—* I i I i 

HP(r), hence (|27[) holds In a similar argument we can prove (|28[) (|29[). 


From (29) we have 
Then 


in the distributional sense. Notice divr'0 and (divr7tcurl M2('0)) are in H 2 (F), then 

□ 


[7iCurlM2(^/’)] =^/’ in H, 2 (F). 


[divr7tCurlM2('0)] = divrt/' 


(30) holds. 
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Now we are ready to prove the equivalence between the transmission eigenvalue problem 


Theorem 2.2 in 15 


and the system of integral equations (24). Our proof follow the lines of the proof of 


Theorem 2.1 The following statements are equivalent: 


(1) There exists non trivial E, Eq G E — Eq G H(curl^,Zl) such that (If)-(ll) 

holds. 


(2) There exists non trivial (M, J)GHj^(r)xH 2 - I (div, T) such that (24) holds and 
either E^(M, J) = 0 or E°°(M, J) = 0 where 


E“(M,J)(h) = hx ( —curl J M{y)e 


-ikx- 


yds,. 


(31) 


1 


iTT^-V J^divr3iy)e-^^^-ydsy + J^Jiy)e 


-ikx- 


yds, 


X X 


with the same expression for E°°(M, J) where k is replaced by ki. 

Proof. Assume (1) holds. Then from the argument above ([2^ we have that M and J 


satisfy (24) and hence it suffices to show E“(M, J) = 0 and E“(M, J) = 0. To this end, 
recall that Eq has the following representation 


Eo(a:) = curl / M{y)^k{x,y)dsy + / 3{y)^k{-,y)dsy 


J dwr3{y)^k{-,y)dsy 


(32) 


where EoXz/ = Exz/ = M and curlEo x u = curlE x o = J. Then, from the jump 


relations (27)-(30) of the vector potentials applied to (32) and (24) (see also [1^), we 


obtain that (Eq x z/)+ = 0, (curlEg x o)+ = 0 (+ denotes the traces from outside of 
D) and hence the far held pattern E“(M, J) varnishes. The asymptotic expression of 


the fundamental solution $(-, ■) in page 23, yields (31). Similarly we can prove that 
E“(M, J) = 0. 

Next assume that (2) holds and dehne 


Eo(x) = curl / M{y)^k{x,y)dsy + / J{y)^ki-,y)dsy 


+ pV J divr3{y)^k{-,y)dsy 


X G 


\r 


with the same expression for E where k is replaced by ki. Again from the jump relations 


of vector potentials and (24) we have 


curl curl E — /c^nE = 0, curl curl Eq — fc^Eo = 0 
E X z/ = Eq X z/, curl E x z/ = curl Eg x z/ 


in D 
on T 
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(note that E and Eq are in L‘^{D). Therefore it suffices to show Eq and E are non trivial. 
Assume to the contrary that Eq = E = 0, and without loss of generality E°°(M, J) = 0, 
then by Rellich’s Lemma (see e.g. 11 ) E = 0 in M.^\D. Hence the jump relations imply 
M = 0 and J = 0 which is a contradiction to the assumptions in (2). This proves the 
theorem. □ 

The above discussion allows us to conclude that in order to prove the discreteness of 
transmission eigenvalues we need to show that the kernel of the operator L(/c) is non¬ 
trivial for at most a discrete set of wave numbers k. 


2.1 Properties of the operator L(A:) 

In the following, we will show the operator L(fc) is Fredholm of index zero and use the 
analytic Fredholm theory to obtain our main theorem. To this end we hrst show that for 
purely complex wave number k := k > 0, L(fc) restricted to 

Hq ^ ~^(div, F) := |u G (div, F), divpu = 0 

satisfies the coercive property. In the following lemma we use the shorthand notation 


Ho(r) ;= H, ^(r)xHo^’ ^(div,r) and its dual space H*(r) := H|(r)x(Ho^’ ^(div,r 


3 __ 

' 2 ’ 2 I 


_3 

where the dual ( Hn ^ ^' 


(div,F)) of the subspace Hq ^ (div, F) C H~S“5(div, F) is un¬ 


derstood in the sense of the duality dehned by Lemma 2.1 


Lemma 2.4 Let k > 0. The operator L(m) : Ho(F) —> H*(F) is strictly coercive, i.e. 


L(m) 


M 

J 


M 

J 


> c 


M 1 -|- J 3 1,1. , I , 

' "H 72 (r) II iiH-3.-2(div,r); ’ 


where c is a constant depending only on k. 


_ _3_ 

Proof. We consider the following problem; for given (M, J) G ^ (F) x Hg ^ (div, F) 
hnd U G L2(M3)^ curlU G L2(M3)^ curCU G such that 


(curl ^ -I- n/t^)(curl ^ k^)U = 0 


m 


^F 


[u X curCU] = {nn^ — on 

[z/ X curCU] = on 


F 

F 


(33) 

(34) 

(35) 


where [■] denotes the jump across F. Multiplying 
by parts yield 


by a test function W and integrating 


/]R3\r 


(curC + nrL^)'\J ■ (curC -|- n^yWdx 


= — K^) 


I 7rcurl W ■ Mds-|- J yrW ■ Jds 


(36) 
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First we show that the right hand side is well dehned. Note that div (curl W) = 0, hence 
from ^ curlW G and thus 7 rCurlW G which implies Jp 7 rCurlW-Mds 

is dehned in ^ (F) duality. Since 7 rW G ^(F) and curlpW = 7 rcurlW G 


2.1 


J^ 7 rW ■ Jds is well dehned. 


(F) then from Lemma 
Now let 

V ;= {U G L2(M3)^curlU G L2(M3)^ curl^U G L2(M3)} 
equipped with the norm 


U||^= / (|curl2U|2 +|curlU|" +lUnda:. 


Next taking W = U in the continuous sesquilinear form in the left-hand side of (36), and 


after integrating by parts (note that U and curl U are continuous across F, we obtain 


(curl ^ -I- nK^)\] ■ (curl ^ -|- K^)\Jdx 


(Icurl Up -I- {nn^ + K'^)|curlUp -I- K^\\]\^)dx > c||U 


V 


where c is a constant depending on k. The Lax-Milgram lemma guaranties the existence 


of a unique solution to (36). Up to here we did not need that divrJ = 0. Next we dehne 


U = curl j - ^^{■,y))ds + j 3{y){^^^{-,y) - ^^{■,y))ds 

+ xn V / divrJ(|/)<l> 7 ;^^(-, 2 /)ds - / divrJ(|/)<F^(-,?/)ds, 


{i^/nKy 


{iny 


then U G L^(M^), curlU G L^(M^), curl^U G L^(]R^) and satishes (33)-(35), hence U 
dehned above is the unique solution to d^. Now for a given yrcurlW G H 2 (F), let us 
construct a lifting function W G H^(M^) 16 such that yrCurlW = yrCurlW, 7 rW = 0 


and ||W||h2(]r 3) < c|| 7 rcurlW|| i for some constant c. Then 
'' ' (r) 


7 rcurl W ■ Mds 


— K^\ 


/R3\r 


7 rcurl W ■ Mds 


(curl ^ -I- ■ (curl ^ Wdx 


< ||U||v||W||v < c||U||v|| 7 rcurl 


Hence ||M|| < c||U||v Similarly for given 7 rW G H 2 (F) we construct the lift¬ 


ing W 2 G H^(M^) 16 such that 7 rW 2 = 7 rW, 7 rcurlW 2 = 0 and ||W2 ||h2(m3) < 
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for some constant c. We recall that divrJ = 0 hence from the Helmoltz 


clK^Woll 3, , 
II ^iiH2(r) 


decomposition J = cnrfrQ' € H 2 (F). Thns we have 


7 rW • Jds 


1 


7 rW 2 ■ Jds 


— ^21 


/R3\r 


(cnrl ^ + nK^)U ■ (cnrl ^ + K^)W 2 dx 


— c||U||v||W 2 IIV < c||U||v||7tW||^3^^^. 

Since J = cmtrh £ (F), then by dnality II J|| 3 1 < c||U||v 

i'i VO j j II iiH-2'-2(div,r) “ " " 

Finally 




M 

J 


M 

J 


I 7 rcnrl U ■ Mds + J 7rU ■ Jds 


> cIlUlIv > c ( ||M 


iv ^ c 1 iiivi|| 1 + J 3 1 ,J. 

IV _ I II Iijj- 2 (r) II llH-2.-2(div,r) 

where c is a constant depending on k. This proves onr lemma. □ 

Next we proceed with the following lemma. 

Lemma 2.5 Let'y{k) := andki = k^/n for k g C\R_. Then L(/i;)+ 7 (/i;)L(i|/c|) : 

Ho(r) —!■ H*(F) is compact. 


Proof. From Theorem 3.8, the operator 

(Sfc, - Sfc) + 7(fc)(S.|fcp - S,|fc|) : H-i(F) ^ Hi(F) 


is compact. Then from (26) we have 


Vr o {Su, - Sk) o divr + 7(fc)Vr o - ^*|fc|) o divr : H -2 (F) ^ H 2 (F) 

{Kk, - Kfc) + 7(fc)(Ki|fcp - K,|fc|) ; H-i(F) ^ H^(F) 

(Kfc, - Kfc) + 7(fc)(K.|fcp - K,|fc|) ; H-5(F) ^ Hi(F) 


Ik,. - Ik,) + -,ik) - ^k.,, ): H-i(r) ^ Ht(r) 


are compact. It remains to show that 

{k.^Sk, - PS,) + 7(fc)((*|fci|)2S,|fcp - (*|fc|)2S,|fc|) : H-i(F) ^ Hi(F) 
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is compact. Since 


- k^Sk) + 7(fc)((z|fci|)2s,|fcp - (*|fc|)2s,|fc|) 

= {k,^{Sk, - So) - k^iSk - So)) + 7(fc)((*|fci|)2(S,|fcp 


(z|fc|)2(S,|fc|-So)) 


and Sfc — So is compact, then the compactness follows. Hence the proof of the lemma is 
completed. □ 


In order to handle the non divergence free part of J, we will split J ;= Q + P where 

3 1 1 I 

Q G Hq ^(div, r), P = VrP G H)*(r) and rewrite the equation (2^ for the unknowns 
(M, Q, P). To this end let us dehne 


Hi(r) := |P G H|(r),curlrP = o} 


and introduce the operator 
L{k) = 


kiTki — kTk 

Kfc, 

-Kfc 



Sfc, 

-Sfc 

Sfcj Sfc 

Kfc, - Kfc 

Sfc, 

-s. 

{Sk,-Sk) + Vro{^Sk, 




o divi 


(37) 


From from Lemma 


2.1 


and Lemma 


2.2 


L{k) : Ho(r) X Hi(r) -G H*(r) X H- 2 (r) is 
bounded and furthermore the family of operators L(fc) depends analytically on fc G C\M_, 


where recall Ho(r) := ^ (F) x Hq ^ (div, F) with its dual H*(F). We hrst notice that 

(24) is equivalent to the following: 

L(fc) 


M 

J 


M 

J 


= 0 , 


for any (M, J) G (F) X H 2’2 (curl, F) which equivalently can be written as 

L{k) = 0, 



for any (M, Q, P) G H* x ^ (F). Now we are ready to prove the following lemma. 

Lemma 2.6 The operator L(A;) : Ho(F) x Hi(F) —)■ H*(F) x H^5(F) is Fredholm with 
index zero, i.e. it can be written as a sum of an invertible operator and a compact operator. 
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Proof. We rewrite the operator L(fc) as follows 

/ l{k){i\ki\Ti\k^\-i\k\Ti\k\) 'y{k)(Ki\k,\ - Ki\k\) 0 

L{k) = - I 7(/c)(Ki|fep - Kj|fe|) 7(fc)(Si|fci| - Si\k\) 0 

\ 0 0 Vr o (-^ + ^)So o divr 

'yi.k') I1 I i|/c|Tj|^|^ 'y{k') 0 

7(^) (Ki|A:i| T(^) (^i|fci| 0 

0 0 


Vr o (-^ + p)So o divr 


— fcTfc — Kfc 

Kfci Kfc Sfc 

Kfci Kfc Sfc 

=: U{k) + U{k) 


(Sfc,-Sfc) + Vro\^Sfc,-^S 


fc2»fc. 


o divi 


(38) 


where Li(fc) is the hrst operator and Li 2 {k) is the sum of the last two operators. Then 
from Lemma 2.5 and the fact that — Sfc,Kfc^ — are s moo thing operators of order 
3,2 respectively, we have Li 2 {k) is compact. From Lemma 2.4 and the fact that So is 


invertible, whence we have Li(A;) is invertible. This proves our lemma. 


□ 


3 The case when N — I changes sign 


In this section we will discuss the Fredholm properties of L(fc) when N is not a constant 
any longer. Our approach to handle the more general case follows exactly the lines of the 
discussion in Section 4 of 


15 , and here for sake of the reader’s convenience we sketch the 


main steps of the analysis. 


3.1 Piecewise homogeneous medium 


To begin with, we assume that D = DiU D 2 such that Di G D and D 2 '■= D \ Di and 
consider the simple case when N = 712 / in D 2 and N = nil in Di where ni > 0, n 2 > 0 
are two positive constants such that (ni — l)(n 2 — 1) <0. Let F = dD, S = dDi which are 
assumed to be smooth surfaces and u denotes the unit normal vector to either F or S 
outward to D and Di respectively (see Figure]^. Let us recall the notations ki = k-^/n{ 
and k 2 = ky/n^. 


For convenience, we let 


s,r 


s r /I K /I 

and ’ be the potentials and given by (18) and (19) 


for densities dehned on S and evaluated on F. The solution of the transmission eigenvalue 
problem (10)-(13) by means of the Stratton-Chu formula can be represented as 


Eo(a;) = curl J {Eq x iy){y)^k{x,y)dsy + J {cnilEo x iy){y)^k{-,y)dsy 
+ / divr(curlEo X i/)(y)<I)fc(-,y)dsy in D 


(39) 
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u 



Figure 1: Configuration of the geometry for two constants 


E(a;) = curl / (E x z/)(y)$fcj (x, + / (curlE x i/)(y)<Ffci(-,y)rfsj, 


J divT(curlE X i/)(y)$fc^(-,y)dsj^ in Di 


(40) 


E(a:) = curl J (E x u){y)^k 2 {x,y)dsy + J {cnrlE x u){y)^k 2 i-,y)dsy 

+ / divr(curlE X i/)(y)<Ffc 2 (-,y)cisy 

^2 Jr 

- curl / (E X u){y)^k2{x,y)dsy - / (curlE x z/)(y)<Ffc 2 (-,y)ds 2 , 


-^V / divT(curlE X i/)(y)<Ffc2(-,y)(iSj^ 

^2 ds 


in Do 


(41) 


Let Exv = EqXv = M, curl Exv = curl Eq x z/ = J on F and Exv = M', curl Exv = i' 
on S. From the jump relations of the boundary integral operators across F and S, we 
have that 


koTl - kTl 


K 2 - K 

K^"—K'" 4-T'" — It'" 

^k2 k2^k2 k^k 


hTl + k,T 


ki 


J_rpE I _L'-pS 
k2 k2' ki ki 


M 

J 

M' 

J' 


Kf/ 

Kl;^ 



J_rpr ,E 
k2 ^2 


(42) 


(43) 


Let us denote by L 2 o(fc), L^’^ik), L 2 i(fc), L^’^(/c) the matrix-valued operators in the 
above two equations in the order from the left to the right from the top to the bottom, 
respectively. By the regularity of the solution of the Maxwell’s equations inside D 2 (see 

), we have (M', J') G ’^(S,div) x ’^(S,div). Then the equation 

M' 

J' 


e.g. 
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L2l(fc) 
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_ 1 _ 1 

where (g, h) G ^(S,div) x ^(S,div) corresponds to the transmission problem 
which is to hnd (E 2 ,Ei) G H;oc(curl,x H(curl,Zi)i) and E 2 such that 


curl curl E 2 — fc 2 E 2 = 0 in 
curl curl Ei — fc^Ei =0 in 
1/ X E 2 — X El = g on 

V X (curlE 2 ) — u X (curl El) = h on 


R\Di 

Di 

S 

E 


and E 2 satisfies the Silver-Mueller radiation condition. By well-posedeness of the trans¬ 


mission problem we have L 2 i(fc) is invertible. Hence pugging in (42) M' and J' from (43) 
we obtain the following equation for M and J 


Hk) 


M 

J 


(44) 


where L(/c) ;= lj 2 o{k) — L^’'"(/c)L 2 i(A:)^ 
we can prove the following theorem. 


^L'"’^(/c). Then in a similar way to Theorem 


2.1 


Theorem 3.1 The following statements are equivalent: 

(1) There exists non trivial E, Eq G L^(Z1),E — Eq G H(curl^,Zl) such that (If)-(ll) 
holds. 


_.±. 3 1 

(2) There exists non trivial (M, J) G ^(T) x H ~2 “ 2 (div, T) such that (44) holds and 
E“(M, J) = 0 where 


E“(M,J)(£) = f X —curl / M{y)e-^'^^'yds 


Ivr 


+ —V / divrJ{y)e-^^^-ydsy + / 3{y)e-^'^^-ydsy ] x 


X 


Now we note that E and T are two disjoint curves and hence we have that L^’^(fc), L'"’^(A;) 
are compact. By writing L(fc) as a 3 x 3 matrix operator L(/c) similar to (38), we can 


have the following lemma directly from Lemma 2.6 


Lemma 3.1 The operator L(/c) : Ho(r) x Hi(r) —>■ H*(r) x H 2(r) is Fredholm with 
index zero, i.e. it can be written as a sum of an invertible operator and a compact operator. 
Furthermore the family of the operators L(fc) depends analytically on k E C\M_. 


This approach can be readily generalized to the case when the medium consists of finitely 
many homogeneous layers. 
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3.2 General inhomogeneous medium 


In a more general case where N = n{x)I in Di, where n G L°°{Di) such that n{x) > a > 0 
but still constant in D 2 , we can prove the same result as in Lemma EH by replacing 
fundamental solution (•,?/) with the free space fundamental G{-,y) of 

AG{-,y) + k'^n{x)G{-,y) = —6y in 


in the distributional sense together with the Sommerfeld radiation condition, where n{x) 
is extended by its constant value in D 2 to the whole space Because ^k 2 (■)!/) ~^{'yy) 
solves the Helmholtz equation with wave number ^2 in the neighborhood of L the mapping 


properties of the integral operators do not change. We refer the reader to Section 4.2 of 15 
for more details. 

In fact the above idea can be applied even in a more general case, provided that N is 
a positive constant not equal to one in a neighborhood of L. More precisely, consider a 
neighborhood O oiV m D (above denoted by D 2 ) with C‘^ smooth boundary (e.g. one 
can take O to be the region in D bounded by V and S := {x — ei/(x), x G L} for some 
e > 0 where v is the outward unit normal vector to L). Assume that N = nl m (9, 
where n 7 ^ 1 is a positive constant, whereas m. D\0 N satisfies the assumptions at the 
beginning of the paper, i.e. is a 3 x 3 matrix-valued function with L°°{D) entries such 
that ^ ■ Re(A^)^ > a > 0 and ^ ■ Im(A^)^ > 0 for every ^ G C^. Then similar result as in 
Theorem 3.1| and Lemma 3.1| holds true in this case. Indeed, without going into details, 
we can express Eq by (39) and E by (41) in O and in H \ (9 we can leave it in the form of 
a partial differential equation with Cauchy data connected to E in O. Hence it is possible 


to obtain an equation of the form (44) where the operator L(fc) is written as 

L{k) = Ln{k) - L^'^{k)A-\k)L^'^{k) 


(45) 


where L„(/c) is the boundary integral operator corresponding to the transmission eigen¬ 
value problem with contrast n — 1 , the compact operators and are dehned 

right below (42) and (43) and A{k) is the invertible solution operator corresponding to 
the well-posed transmission problem 


curl curl E2 — k^n2Ei2 = 0 
curl curl El — k'^N'Ei = 0 
i/xE 2 — i^xEi = g 
V X (curlE2) — u X (curl El) = h 


in 

R\Di 

in 


on 

E 

on 

S 


and E 2 satishes the Silver-Miiller radiation condition. 


to prove analogues Theorem 3.1 and Lemma 3.1 


Hence the above analysis can apply 


For later use in the following we formally state the assumptions on N (here (9 is a 
neighborhood of T as explained above). 


Assumption 3.1 N is a 3 x 3 symmetric matrix-valued function with L°°{D) entries 
such that f ■ Re{N)f > a > 0 and f ■ Im{N)f > 0 for every f G |.^| = 1 and N = nl 
in O where n ^ 1 is a positive constant. 
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4 The existence of non transmission eigenvalue wave 
numbers 


In this section we assume that N satisfies Assumption 3.1 and consider pure imaginary 
wave numbers k and, for convenience, let A := —/c^ be a real positive number and start 
by proving an a priori estimate following the idea of (2^ for the scalar case. 


Lemma 4.1 Assume that N satisfies 3.1 and x(a;) G C“(iA) is real valued cutoff function 
with 0 < X < 1 and y = 1 m D\0. If'v^ and 

(curl curl + A)v = 0 in D 

then there exists a constant K (y) such that for sufficiently large A 

11(1-X)vf 


A 


llxviP < K 
Proof. Since y G C“(iA) we have 

0 = / (curl curl + A)v ■ {x^v)dx = / curl curl v ■ {x‘^v)dx + A / v ■ {x‘^v)dx 

J D Jd J d 


(46) 


curl V • curl (x^v)(ix + A / v ■ {x^'v)dx 


' D 


>D 


/ curl V ■ (ycurl (xv))(ia; + / curl v ■ (Vy x (xv))(ix + A / v ■ (y^v)(ix 
Id Jd Jd 

/ curl (yv) ■ curl (yv)(ia; — / curl (yv) ■ (Vy x \')dx 
Id Jd 

/ curlv ■ (Vy x (yv))(ix + -^ / ^ ■ {x^v)dx 
Id j d 

/ |curl (yv)l^(ia; — / (ycurl v + Vy x v) ■ (Vy x v)(ix 


>D 


>D 


/ curlv ■ (Vy x (yv))(ia; + A / v ■ {x^\')dx 
Id Jd 


I curl (yv)p(ix 


|(Vy X v)p(ia: + A / lyv^da: 


' D 


'D 


' D 


((ycurl v) ■ (Vy x v) — (ycurlv) ■ (Vy x v)) dx. 


>D 


Taking the real part yields 


/ |curl (yv)l^(ix + A / \x'v\‘^dx = / |(Vy x v)|^(ia; 
Id Jd Jd 


and then 


A||yv||2 < A:(y)||vf < K{x) (||xv||^ + 11(1 - x)v|P) 
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which yields (46) for sufficiently large A. 


□ 


Now we are ready to prove the following theorem. 


Theorem 4.1 Under the assumption \3.1\ there exists a sufficiently large real A > 0 where 
A = —k"^ such that (10)-(13) has only trivial solutions. 


Proof. Assume first n — 1 < 0 in C>, let u = E — Eq G Ho(curl^, D), v = AEq G L^(iA), 
then 


curl curl u + AA^u =—(A^ —/)v in D 

curl curl v + Av = 0 in D 

u X u = u X (curl u) = 0 on F. 


Then for any cp G C“(Z1), interpreting (48) in the distributional sense yields 


(47) 

(48) 

(49) 


v(curl curl ip + Xip) = 0, 


>D 


and hence the denseness of C“(Z1) in Ho(curl ,Z1) (see 16 ) yields 


/ V ■ curl u + A / v-u = 0 
'D J D 


Multiplying (47) by v yields 


(50) 


/ V ■ curl + A / A^u ■ \dx + {N — J)v ■ vdx = 0 
Id Jd Jd 


Combining the above with (50) yields 


X {N — I)u ■ \dx + {N — /)v ■ 'vdx = 0 
Jd J d 


Multiplying (47) by u and integrating by parts yields 


(51) 


/ Icurlupda; + A / A^u ■ uda; + / (A^ —/)v ■ udx = 0 
Jd j d Jd 

Noting that N is symmetric, we have (A^ — /)u ■ v = (A^ — /)v ■ u and hence 


/ |curl u\‘^dx + X Nu ■ udx + {N — I)u ■ \dx = 0 
Id Jd Jd 


(52) 


By regularity 26 v is sufficiently smooth in D away from the boundary and hence by 
unique continuation we can see /(^(n — 1)(1 — x^)|vpda; 7 ^ 0. Then combining (51) with 
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(52) yields 


>D 


|curl + A / Nu-udx=— j {N — I)^ ■ v dx 


1 

A 


' D 
. 2 ,, — 


' D 


(iV —/)x^v ■ V dx + / (iV —/)(1 — x^)v ■ vdx 


' D 


' D 


If, ,, 9s _ f fn(^ — • vdx 

A./z)' A y ^ fjN-I)(l-x^)v-vdx 


= vin - 1) (1 - X^)H^ dx 1 + 


A 


'O 


fpi^ - -vdx 

(n-l)fa(^-X^)H^dx 


From Lemma 4.1 we have for sufficiently large A 

vdx| K(Nmax + l) 


(l-^)/o(l-A^)|v|2dx 


< 


< 1 , 


where N^ax is supremum over D of the largest eigenvalue of N, which implies 

/ fn(-^ ~ d)x^V ■ vdx , 

»(1 + , •'°A f2'| |„ I > 0. 


(53) 


(54) 


Then, since n — 1 < 0, the real part of (54) is non positive for sufficiently large A but 


the real part of (53) is non negative. Hence the only possibility is u = 0,v = 0, i.e. 
E = Eo = 0. 

Let us next consider n — 1 > 0 in (9, and let u = E — Eq, v = AE. Then 


S V 

= AE. Then 


in 

D 

(55) 

in 

D 

(56) 

on 

T 

(57) 


Using the same argument as for (50) 


/ curl ■ vdx + A / iVv ■ udx = 0 
Id J d 


Multiplying (55) by v yields 


(58) 


/ V ■ curl ^udx + A / v ■ udx + (N — /)v ■ vdx = 0 
Id Jd Jd 


Combining the conjugate of the above with (58) yields 


X {N — J)u ■ vdx = / N — Iv ■ vdx 
Jd Jd 


( 59 ) 
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Multiplying (55) by u and integrating by parts yields 


|curlu|^(ix + A / |u|^(ia:+ / {N — ■ udx = 0. 


>D 


>D 


' D 


Note that since N is symmetric, then (A^ — /)u ■ v = {N — J)v • u and hence 


/ Icurlul^dx + A / \u\'^dx+ / (A^ — J)u ■ vda; = 0. 
Id Jd Jd 


Then combining (59) with (60) yields 


( 60 ) 


/ |curlu| dx + A / |u| da; = — - / N — I v ■ v dx 
Id Jd ^ Jd 

1 


1 I y A^ —/v-vda;+ / (1 — y )A^ — Jv ■ vda; 

\ Jd Jd 

Id — / V ■ vdx 


= —— / (1 — y )A^ — Jv • vdx 1 + 


A 


' D 


= / (n-l)(l-y2)|v|2dx 1 + 


A 


'o 


f^(l — y2)A^ — /V ■ vdx 
Id — Iw - w dx 


(61) 


(62) 


From Lemma 4.1 we have for sufficiently large A 


|/By^^-^v-vda:| ^ K{N^a^ ^ c 1 


- X^)W? dx 


A 


Then 


I fn x'^dX — Iw -w dx 

^ I ^ + ~T~( -TATt -2p 12^ I ^ 

J^(n- 1)(1 -y^)|v|^da: 


Therefore, since n — 1 > 0, the real part of (62) is non positive for sufficiently large A 
but the real part of (61) is non negative. Hence the only possibility is u = 0,v = 0, i.e. 
E = Eq = 0. □ 


5 Discreteness of transmission eigenvalues 

Recall that in Section]^ we have proved that L(/c) is a Fredholm operator. Hence to show 
discreteness we will use the analytic Fredholm theory 11 . To this end we must show that 
there exists k such that L(A;) is injective. 


Lemma 5.1 Assume that N satisfies 3.1. There exists a purely imaginary k with suffi¬ 


ciently large \k\ > 0 such that L(A;) is injective. 
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Proof. Let us extend N to M.^\D hj N = nl where n is the constant N\o- Assume there 
^ such that L(/c) ( ^ ] =0. We will show that if k is purely imaginary 


exists 


with large modulus, then 


M 

J 


= 0. Recalling (45), we dehne 


M' 

J' 


= A ^{k)L^’^{k) 


M 

J 


and 


Eo(a;) = curl J M{y)^k{x,y)dsy + J 3{y)^ki-,y)dsy 

+ J divT3iy)^k{-,y)dSy in R^\T. 

From the dehnition of ( ^ ] there exists E G L(Zi)i), Di := D\ O, such that 


Also we dehne 


curl curl E — k'^N'E = 0 in Di 

[E X z/]+ = M' on S 

[curlE X //]■*■ = J' on S. 


E(a;) = curl / M{y)^k 2 {x, y)dsy 


Jiy)^k2{-,y)dsy 


/ divTJ(y)<hfc2(-,y)dsy 

^2 Jr 

curl / M'{y)<^k 2 {x,y)dsy - / J'(y)<Ffc 2 (-,y)dsj^ 


■^y7 j divT3'iy)^k2i-,y)dSy in 




Jump relations across F applied to E, Eq along with the equation (44) yield 

curl curl E — k'^N'E = 0 


curl curl Eq — /c^Eq = 0 

{v X E)^ = (i/ X Eo)^ 
(z/ X curlE)^ = (z/ X curlEo)^ 


m 

in 

on 

on 


^3\F 

^3\F 


F 

F. 


(63) 

(64) 

(65) 

( 66 ) 


From Theorem 4.1 if fc is purely imaginary with large enough modulus then ([6^ -([6^ in 
D only has trivial solutions. Since N = nl where n is a constant in M.^\D, then the 
variational formulation of (63)-(66) in M.^\D is (36) where the right hand is 0 and ]R^\F 
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is replaced by M.^\D. Then U 
(|2^-® yield M = 0 and J = 


= 0 and hence E = 0, Eq = 0 in ]R^\r. The jump relations 
0 and this proves the lemma. □ 


Finally, combining Lemma 3.1 and Lemma 5.1 


orem using the analytic Fredholm theory 11 . 


we can immediately prove our main the- 


Theorem 5.1 Assume that N satisfies Assumption 3.1 
eigenvalues in C is discrete. 


Then the set of the transmission 
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